Ann. Phys. (Leipzig) (1999) ,1-6 



Metallic phase between the Fermi glass 
and the Wigner crystal in two dimensions 

X. Waintal 1 , G. Benenti 1 , and J.-L. Pichard 1 

1 CEA, Service de Physique de l'Etat Condense, 

Centre d'Etudes de Saclay, F-91191 Gif-sur-Yvette, France 

pichard@spec . saclay . cea . f r 

Received 21 July 1999, accepted dd.mm.yyyy by ue 

Abstract. For intermediate Coulomb energy to Fermi energy ratios r 3 , spinless fermions 
in a two-dimensional random potential form a new quantum phase, different from the Fermi 
glass (weakly interacting Anderson localized states) and the Wigner crystal (regular array of 
charges pinned by the disorder). The intermediate phase is characterized by an ordered flow 
of persistent currents with a typical value decreasing with excitation energy. Extending finite 
size scaling analysis to the many body ground state, we find that electron-electron interactions 
can drive the Fermi glass towards an intermediate metallic phase (Coulomb metal). 
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1 Introduction 

According to the scaling theory of localization jl|], all states are localized for non- 
interacting particles in two dimensions (2d). This absence of a metallic state in 2d is 
nowadays challenged since a lot of transport measurements, following the pioneering 
experiments of Kravchenko and co-workers Q , give evidence of an insulator- metal tran- 
sition (IMT) when the carrier density is increased. Many experimental results suggest 
that Coulomb repulsion drives this phenomenon: (1) In a very clean heterostructure 
the IMT was observed at a Coulomb energy to Fermi energy ratio r s ss 35 ||], close 
to r s i=s 37 for which a clean Wigner crystal melts according to Monte Carlo simula- 
tions (2) In more disordered samples, the IMT typically occurs at r s ss 10, again 
compatible with the value expected for the melting of a pinned Wigner crystal (3) 
At a weaker r s interactions no longer dominate and a re-entry towards an insulating 
phase (Fermi glass of weakly interacting localized particles) was observed at r s ss 6 in 
||. As at the IMT kpl sa 1 (hp denoting the Fermi wave vector and I the elastic mean 
free path), disorder is important and has to be considered, together with interactions, 
in a non perturbative way. 

Taking advantage of numerical techniques, we consider a model of N Coulomb 
interacting spinless fermions in a disordered square lattice with L 2 sites. The Hamil- 
tonian reads: 

H = - t E c i^+E^ + c/ E^ (!) 
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where c\ (cj) creates (destroys) an electron in the site i, t is the strength of the 
hopping terms between nearest neighbours (t = 1 in the following) and the inter- 
particle distance for a 2d torus. The random potential v% of the site i = {i x , i y ) with 
occupation number n, = cjcj is taken from a box distribution of width W. The 
interaction strength U yields r s — U / (2ty/im e ) for a filling factor n e — N/L 2 . A 
Fermi golden rule approximation for the elastic scattering time leads, for n e < 1, to 
k F l w 192Tm e (t/W) 2 . 



2 Persistent currents 

To measure derealization effects induced by interactions, we study the sensitivity of 
the ground state to a change in the boundary conditions . Boundary conditions are 
always taken periodic in the transverse y-direction, and such that the system encloses 
an Aharonov-Bohm flux 4> in the longitudinal x-direction. Imposing <f) — n/2 ((f> = it 
corresponds to anti-periodic boundary conditions), one drives a persistent current of 
total longitudinal and transverse components given by 

with E ground state energy. The local current l\ flowing at the site i in the lon- 
gitudinal direction is defined by l\ — 2Im(>I , ol c I x +i i s c i x % v l^o)j with an analogous 
expression for If (l^o) is the ground state wavefunction). The response is paramag- 
netic if Ii > and diamagnetic if /; < 0. Exact diagonalization techniques (Lanczos 
method) are possible only for small system sizes and here we consider = 4 particles 
in L 2 = 36 sites, for W = 5, 10, 15 (k F l = 2.7,0.7,0.3 respectively) and < r s < 42. 
We summarize in Fig. results from a statistical study of an ensemble of 10 3 clusters. 
\Ii\ and \I t \ have acceptable log-normal distributions for all values of r s when W > 5 
[0 and therefore the log-averages shown in Fig. |l| can be considered as typical values 
of the persistent currents: I t = exp < ln|Jj| >, with an analogous expression for 
-fjtyp (brackets indicate ensemble average). Fig. [l] (a) shows that the ground state 
transverse current is suppressed before the longitudinal current. It is then possible to 
single out an intermediate phase, in which the local currents flow in an ordered way 
along the direction imposed by the flux [Q, |). This behaviour is evident in Fig. [j] (b): 
the ground state local current angles Qi = arctan(7|/J,-) undergo a transition from a 
nearly uniform angular distribution (glass of currents randomly directed in the plane) 
towards a distribution strongly peaked in the longitudinal direction. As a consequence, 
the sign of the magnetic response becomes independent of the microscopic realization 
of the random potential (see Fig. |l| (b) insert). Notice that, for a not too strong dis- 
order (hp I > 1), in the intermediate phase J; is enhanced by interactions up to a 
factor 3 (for W = 5, r s ps 6) before being suppressed due to charge crystallization 0. 
Fig. H (c) displays the behaviour of the longitudinal current for the 20 lowest energy 
levels. The enhancement in 7 ; ^yp disappears away from the ground state. In addi- 
tion, it is possible to identify three different regimes when W = 5: for approximately 
r s < 2 and r s > 17 the longitudinal current increases with the excitation energy e, all 
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Fig. 1 (a) Log-averages of the ground state longitudinal current 7; (open circles) and 
transverse current It (filled circles) at W = 5. Insert: the same at W — 15. (b) Distribution 
Pe of the ground state local currents angles at W = 5, for r s = (dotted line), r a = 6.3 
(full line), and r s = 42 (dashed line). Insert: fraction c p of paramagnetic samples for W = 5 
(circles), W = 10 (squares) and W = 15 (diamonds), (c) Log-average of the longitudinal 
current at W = 5, for the ground state (circles), the first excited state (squares), the second- 
fourth excited states (diamonds), the fifth-ninth excited states (triangles up), and the tenth- 
nineteenth excited states (triangles down). Arrows indicate increase (or decrease) of the 
longitudinal current with excitation energy. Insert: the same at W = 15. (d) Same as in 
(c) but for the transverse current at W = 5. Insert: fraction of paramagnetic samples as a 
function of the excitation energy e at W = 5, for r s — (circles), r s = 6.3 (squares), and 
r s = 42 (diamonds). 
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the contrary for intermediate r s values. We remind that experimentally the metallic 
(insulating) phase has been identified (2[ [|, |(| from a decrease (increase) with tem- 
perature of a related transport property, the conductance. For W = 15, / ; ^yp always 
decreases with e (Fig. ^ (c) insert), indicating that a too strong disorder probably 
forbids a metallic phase in the thermodynamic limit. Fig. [j] (d) shows that, increasing 
e, the transverse current is suppressed at larger r s values, eventually compatible with 
r s values at which the longitudinal current is suppressed. This fact suggests that the 
intermediate phase is destroyed at large enough temperatures. For e > 2 (the Fermi 
energy is given by 6f — %n e t — 0.35) the sign of the persistent current is disorder 
dependent at any r s value (Fig. Ill (d) insert). 



3 Finite size scaling for the N body localization length 

Though small clusters exhibit an enhancement of for intermediate r s , exact diago- 
nalization does not allow to vary system size and to establish if the intermediate phase 
is metallic at the thermodynamic limit. In order to answer this question and to un- 
derstand how Coulomb interaction destroys Anderson localization, we extend finite 
size scaling analysis [|l0| [llj to interacting systems. We consider an ensemble of 5 x 10 3 
clusters with iV = 3, 4, 5 particles in square lattices of size L = 24, 28, 31 respectively, 
corresponding to very low filling factors n e ss 5 x 10~ 3 . To have Anderson localization 
inside these sizes we consider a large disorder to hopping ratio W/t — 10. Therefore 
the low energy tail of the one body (IB) spectrum is made of impurity states trapped 
at some site i of exceptionally low vi. As we are interested in studying the effect of 
Coulomb repulsion on genuine Anderson localized states we get rid of the band tail. 
Typically we ignore the L 2 /2 first IB levels (but results do not change provided that 
the Fermi level if out of the band tail, ep > — At; under this condition one can roughly 
estimate kpl « 1). From this restricted subset of IB states we build a basis for the N 
body (NB) problem, truncated to the Nh = 10 3 Slater determinants of lowest energy 
(convergence tests are discussed in || ) . The scaling ansatz for the NB problem reads 




where we assume it is possible to map the localization length £j, at the system size L 
onto a scaling curve f(L/£), with £ characteristic scaling length of the infinite system. 
To characterize the NB ground state by a suitable localization length £l , we consider 
the change 5pj of the charge density induced by a small change Svi of the random 
potential Vi located at a distance \i — j\. To improve the statistical convergence, we 
calculate more precisely the change Sp(r) = J2j ^Pr,j v of the charge density on the L 
sites of coordinate j x — r yielded by the change t?o,i v — > l-01uo,i„ for the L random 
potentials of coordinate i x = 0. Fig. ^ (a) shows that |<5p(r)| is reasonably fitted by a 
log-normal distribution. Therefore it makes sense to characterize the typical strength 
of the fluctuations by £>/Otyp( r ) = CX P < ni l^/ 9 ( r )l > an d extract the length £l over 
which the perturbation is effective from the exponential decay <5/9typ( r ) K ex P(~ r /£-L)- 
Such a decay occurs only over a scale r << L/2 since the boundary conditions are 
periodic. In Fig. || (b) the lengths £l are obtained from the slope of the linear parts 
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Fig. 2 (a) Distribution of the logarithms of the charge density variation \Sp\ for L — 31, 
r a — 0, W = 10, r — 5, fitted by a log-normal function of mean p = —8.6 and variance 
a 2 — 4.2. Insert: the same at r s = 5.9, with p = —7.5, a 2 = 3.0. (b) Log-average of \Sp\ as a 
function of r, for L — 28, W = 10, r s = (circles), r s = 2.0 (squares), r s = 3.9 (diamonds), 
and r s = 5.9 (triangles). The lengths £z, are given by the inverse slopes of the straight lines. 
Insert: the same for a quasi-one dimensional geometry (L x — 80, L y = 8) at W — 8 (ensemble 
of 10 3 clusters), (c) Reduced localization length £l/L as a function of r B for 14 7 = 10, L — 24 
(circles), L = 28 (squares), L = 31 (diamonds). Insert: ratios £l/L mapped onto the scaling 
curve / as a function of the reduced scaling length Localized branch (filled symbols): 

r s = (triangles up), r 3 = 1.0 (triangles left), r s = 2.0 (triangles down), r 3 = 2.9 (triangles 
right), r s — 3.3 (circles), r s = 3.9 (squares). Delocalized branch (open symbols): r 3 = 4.5 
(squares), r s — 4.9 (triangles), r s = 5.3 (diamonds), r s = 5.9 (circles), (d) Log-average of 
\5p\ as a function of r s for N = 4 particles in a 6 x 6 lattice, r = 1, W = 5 (circles) and 



W = 15 (squares) (ensemble of 5 x 
lowest energy one. 



10 clusters). IB states are filled in starting from the 
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of the curves (straight lines). A good exponential decay on a much larger interval can 
be observed in a quasi one dimensional geometry (L x 3> L y ). Fig. || (c) gives how 
the reduced localization length £l/L depends on r s for the three considered sizes. A 
critical point appears at rf w 4.3: for r s < rf , £l/L decreases with L (Fermi glass) 
while it increases for r s > rf (Coulomb metal). Fig. || (c) insert verifies the scaling 
ansatz (||): assuming suitable scaling lengths £ all the data can be mapped onto a 
single scaling curve /, which develops two branches, as typical of second order phase 
transitions (like the IB Anderson transition in three dimensions jlO[ |llj). A power fit 
of £ oc \r s — rf \~ u yields a rough estimate for the critical exponent v w 4. In the quasi 
one dimensional case we were able to detect a strong interaction induced enhancement 
of the localization length (see Fig. || (b) insert) but no signature of a phase transition. 
Fig. |^ (d) shows that a local perturbation in the random potential plays a role similar 
to a change in the boundary conditions. For the small system sizes accessible to exact 
diagonalization it is not possible to extract the length However, it is possible to 
see how the typical value <5/%yp( r )j evaluated for example when r = 1, evolves with 
r s . This quantity increases in the intermediate phase before decreasing due to charge 
crystallization. Unfortunately it is not possible, within the approximate technique 
described above, to evaluate how £l changes with L in the Wigncr crystal phase. 



4 Conclusions 

Coulomb repulsion can drive a two dimensional system of spinless fermions in a random 
potential towards a new metallic phase, the Coulomb metal, which is different from 
the Fermi glass and the Wigner crystal. The local current could be the order 
parameter driving the Fermi glass-Coulomb metal transition (through Oi) and the 
Coulomb metal- Wigner crystal transition (through 
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